Coulomb scattering in a 2D interacting electron gas and production of EPR pairs 
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We propose a setup to generate non-local spin-EPR pairs via pair collisions in a 2D interacting 
electron gas, based on constructive two-particle interference in the spin singlet channel at the n/2 
scattering angle. We calculate the scattering amplitude via the Bethe-Salpeter equation in the ladder 
approximation and small r s limit, and find that the Fermi sea leads to a substantial renormalization 
of the bare scattering process. From the scattering length we estimate the current of spin-entangled 
electrons and show that it is within experimental reach. 

PACS numbers : 73.23.-b, 71.10.Ca, 3.67.Mn, 



In recent years, the spin degree of freedom of electrons 
has become of central interest in semiconductor research 
PJ. This is particularly so for spin-based quantum in- 
formation processing, where the basic resources are EPR 
pairs such as non-local spin singlets formed by two elec- 
trons that are spatially separated 0. A number of recent 
publications have described ways to produce such spin- 
correlated two-electron states, as well as orbital entangle- 
ment 0. Here, we propose a new concept based 
on a two-particle interference mechanism that is well- 
known from elementary scattering theory Q: the cross- 
section for two electrons in vacuum, given in terms of the 
scattering amplitude / by \ S / T {9) — \f(0) ± f(n — 9)\ 2 , 
favors singlet (+) over triplet (— ) states in the out- 
going channel around the scattering angle — ir/2. 
Thus, in principle, two-particle scattering processes can 
be used to generate EPR pairs. However, in the context 
of solid state systems, the question immediately arises 
whether this scattering effect remains operational, and 
moreover observable, in the presence of a Fermi sea con- 
sisting of many interacting electrons such as a typical 
two-dimensional electron gas (2DEG) formed in GaAs 
heterostructures. In this Letter, we will show that within 
Fermi liquid theory the answer is affirmative. 

We focus on 2DEGs since these systems are promising 
candidates for the observation of such effects. Indeed, re- 
cent experiments 0] have demonstrated that in a 2DEG 
the flow of electrons as well as scattering off impurities 
can be controlled and monitored via AFM technology. 
Thus, we believe that a setup as shown in Fig. ^a) is ex- 
perimentally realizable and should allow the observation 
of the angular and density dependence of the scattering 
cross section. Once the EPRnairs are created their sin- 
glet character can be tested [S( by a noise measurement 
in a beam s plit ter configuration f2j, or by tests of Bell 
inequalities |lf| . 

We note that the experimental observation of this type 
of entanglement would provide support for the appli- 
cability of Fermi liquid theory to 2D systems, which 
has been questioned by Anderson fTl) . Indeed, one 




Figure 1: (a) Proposed setup: two quantum point contacts 
(QPC) filter electrons from two reservoirs with initial mo- 
menta pi ~ — P2. The two detectors (with an aperture angle 
259) are placed such that only electrons which collide (shaded 
area) at a scattering angle around 7r/2 are registered. Because 
of interference, the scattering amplitude / vanishes at tt/2 
for the spin-triplet states, allowing only the spin-entangled 
singlets to be collected: one electron of the singlet state in 
detector 1 and its partner in detector 2. The scattering cross- 
section and the electron flux could be measured via an AFM 
tip ,8J. (b) Scattering parameters: P = pi + P2 = p'i + P2 is 
the total momentum, p = (pi — P2)/2 and p' = (p'i — P2)/2 
are the relative momenta, and 9 = Z(p, p') is the scattering 
angle between them. The initial (pi,P2 ) and final (p'i,P2) 
momenta are connected by a circle of radius p = p' due to 
energy and momentum conservation. 

can hardly imagine singlet pairs of particles that are 
separated by mesoscopic distances without well-defined 
fermionic quasi-particles. 

For electrons incident from unpolarized sources, we ex- 
pect the ratio of singlets \S) to triplets \Ti), i = 0, ± to 
be 1:3, this mixed state being described by 1/4 IS) (S\ + 
l/4^ i=0 ± \Ti) (Ti\. Our goal is to calculate the ratio, 
of scattered triplets to scattered singlets, for an aperture 
angle 280 of the detectors around 9 — n/2. For small 
<56»this ratio is R ~ S9 2 \f/f\ s=7r/2 , where / is the many- 
body scattering amplitude and /' its angular derivative. 
Solving the Bethe-Salpeter equation in the small r s limit 
we find that |/'//| 2 ~ 1 at = ir/2. This means that 
for small SO the scattering process produces dominantly 
singlets in the direction of the detectors 1 and 2; see 
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Fig. nja). As an experimental check, this ratio can be 
increased by reducing the amount of singlets in the in- 
coming channels, which in turn can be achieved e.g. by 
using spin polarized electron sources or devices that act 
as spin filters such as a quantum dot [3] or a quantum 
point contact (T3 |. 

Setup. We begin with the description of the setup, 
shown in Fig. HJa). Electrons escaping from thermal 
reservoirs with momenta pi ~ — P2 are filtered by two 
quantum point contacts (QPC) and injected into a 2DEG 
where they scatter off each other. To collect after the 
collision only the entangled singlets (EPR pairs), we 
place two detectors such that only collisions of electrons 
with final momenta pi,p 2 and with scattering angle 
6 G [n/2 — 86, n/2 + 58] are registered. Below we estimate 
the expected singlet flux (current). 

As seen from Fig. Qlb), the energies of both incident 
electrons need to be known, in general, to determine the 
scattering angle 6 ~ Z(p',p). However, for the special 
case with opposite momenta p2 — —pi, is easily deter- 
mined by 9 ~ Z(pi,pi). Moreover, the energies are then 
individually conserved, pi ~ p 2 — Pi — P2 (Pi = IPj'Dj 
which ensures that the outgoing scattering states are 
unoccupied (p' 12 > &f)- Finally, to have well-defined 
quasiparticle states with long lifetimes, we assume the 
electrons to be injected with small excitation energies 
& = h 2 p 2 /2m - E F < E F , where E F = h 2 k 2 F /2m is the 
Fermi energy and m the effective mass. 

Scattering t-matrix. Let us evaluate the scattering t- 
matrix for two electrons in the presence of the Fermi sea. 
The condition p 2 = — pi defines the Cooper channel, and 
thus we can follow the work by Kohn and Luttinger on 
interaction-induced superconductivity in a 3D Fermi liq- 
uid In contrast to Ref. ^ij], we consider here a 
2D system where the screened Coulomb potential is non- 
analytic. Moreover, we need the complete angular de- 
pendence of the scattering amplitude — rather than only 
the asymptotics of its Legendre or Fourier coefficients. 
The i-matrix can be obtained from the (direct 1 ) vertex T 
governed by the Bethe-Salpeter equation 0, B3| 



The bare Coulomb interaction in 2D is given by 



T{p',p;P) = A(p',p;P) 

+ h(ky J d ~ k Hk,p;P)G(ki)G(k2)r(p',kP) 



(i) 



where G is the single-particle Green function and A the 
irreducible vertex (we suppress the spin indices as spin is 
conserved). The relative and total 4- momenta are de- 
noted by p = (pi — p%)/2, and P = p\ + P2, where 
Pi = (Pii^i) with the frequencies Wj. The full vertex for 
singlet /triplets contains the direct and exchange parts, 
i.e. T(p',p;P)±T(-p',p;P), 
tain the i-matrix: t — T(u>i ^ 
from it the (unsymmetrized) scattering length A = I f I 2 
via the 2D scattering amplitude / = —tm/Ti 2 ^2irp |lfi| . 
The scattering lengths for singlet /triplets are then given 
by \ S /t(6) = \f(6)±f(e-n)\ 2 . 



From the vertex we ob- 
wa & + 6) 0, and 



^c(q) = ^ e o/l with e, 



'/4neoe r and the dielectric 



constant e r . In a first stage, we use the RPA approx- 
imation to account for screening by the Fermi sea 
this yields a renormalized G and a screened interaction 
V(q) = Vc(<3 , )/[1 — Vc(q)x (9)] given in terms of the bub- 
ble susceptibility diagram \° with q = fa, ui q ) = p'—p the 
momentum transfer. We can consider [13] the usual static 
Thomas-Fermi screening V(q) = e\j{q + k s ), with the 
screening momentum k s — 2me 2 ) /Ti 2 and r s = k s /kF\/2. 
Note that RPA requires a high density, r s <C 1. 

In addition to the single interaction line V(q), the ir- 
reducible vertex A contains, in lowest order in V, two 
diagrams: the crossed diagram and the "wave function 
modification" shown in Fig. 2(c) and (d) of Ref. 0|. It 
is easy to see that these are smaller than V by a factor 
r s jig}. Altogether, this justifies the ladder approxima- 
tion, which consists in keeping only the single interac- 
tion line in A ~ V. Note that the ladder approxima- 
tion requires Xk F -C 1 H^- This condition is consistent 
with RPA: e.g., in the Born approximation t ~ Vc 
Xk F ~ [Vc(k F )m] 2 /2irh 4 ~ (k s /k F ) 2 . Finally, we can 
approximate G by the free propagator Go [l4| . We start 
with the zero temperature T — case, and discuss finite 
T effects later. 

Before solving Eq.Q to all orders, we first consider 
its second-order iteration r( 2 )(p', p; P) = V(p' — p) + 
i/h{2nf Jd~kV(k - p)G(k 1 )G(k 2 )V(p' - k). The lu- 
integration of the Green functions yields 
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= :D(ki,k 2 ), (2) 



with 7V(ki, k 2 ) = 1 - n{k\) - n(k 2 ), n(k) = 6(-^), and 
the P-frequency is HQ — * ^1 + ^ 2 to obtain the i-matrix. 
We now take p = Pi = P2, k = ki = — k 2 so that q — 
|p' — p| = 2k F \ sin 0/2 1 . This yields a single sharp edge 
in the numerator when = 0, which coincides with the 
zero of the denominator at £fe = £ for incident electrons 
with vanishing excitation energies £ = £i. 2 — > 0. Thus, 
the main contribution to the energy integration comes 
from virtual states at the Fermi surface, i.e. ~ 0. We 
set k = k F in V, and integrate only on D(k) = D(k, k). 
The dominant contribution comes from the vicinity of 
the Fermi sphere k — k F . Introducing the notation v := 
(1/2tt) J™dkkD(k), we obtain 



m £, 
— log — 



2ttH 



E F 



(3) 



For finite temperatures with £ <C fcsT 1 *C E F , the oc- 
cupation function n(k) = (1 + e 5fc /' CBT )~ 1 cuts the log 
divergence and £ is replaced by ksT in Eq. 10. This 
logarithmic divergence reveals a 2D Cooper singular- 
ity (see below) very much like in 3D Next, we 



examine the situation slightly away from the Cooper 
channel, i.e. for non-vanishing total momentum. In- 
deed, for experimental reasons it is preferable to have 
a small but finite angle 2a — Z(pi,— P2) between 
the incident particles to avoid flux misalignment with 
no collision at all. In this case, at T = and for 
Pi = P2, we find that N(k, </>)~8(fc — kp — pa\ sin^|) — 
Q(kp — k — pa\ sin 0|) depends on the integration an- 
gle 4> — ^(k, p) and, as a consequence, we find now 
h>(tfi) ~ (m/2-Kh 2 ) log (2a\ sin0|) in the limit p — » k F . 
Thus the cut-off of the log-divergence is determined by 
max{a| sin(j)\, ^/ Ep, ksT/Ep}. However, we recall that 
both outgoing scattering states must be unoccupied, 
which at T = requires p' 12 > &f- For 7r/2-scattering, 
Pi 2 = Pi (cos a ± sin a) p\ > kp(l + a); see Fig. ^b). 
Therefore, a < £,/Ep and the relevant cut-off is still given 
by £ or k B T . 

We repeat this procedure in each order in V and 
rewrite Eq. (JTJ as t(0)=v(0)+(v/2ii) f d<jru(<j>)t(0 - 0), 
with v((j})=2ne 2 /(2k F | sin </>/2| + fc s ) the 2D Coulomb po- 
tential V at the Fermi surface. To solve this equation, we 
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expand v and t into Fourier series: v((f>)=J2 n v ^ e 
We finally get for the i-matrix at the Fermi surface 



t{0) 



^ 1 



with the Fourier coefficients 
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(5) 



and 7 = arcsin (r s /V2). Below we use this result to 
evaluate the scattering length. 

Scattering length. To illustrate how the scat- 
tering process gets renormalized by the Fermi sea, 
we compare the above result Q with the i-matrix 

|tc(0)| = e ?i[tanh( /k F h 2 ) Tr/mkp sin 2 (0/2)} ' 
obtained for the bare Coulomb potential Vc in 2D 
We use typical parameters for a GaAs 2DEG: e r = 13.1, 
r s = 0.86 and a sheet density n — 4-10 15 m -2 and as- 
sume £ < k B T = IQ- 2 E F (T = 2 K). In Fig. EJa) we plot 
the scattering length X(0) = \tm / h 2 ^/2ttp\ 2 as function 
of the scattering angle (without antisymmetrization) . 
The reduction in amplitude due to the Fermi sea is seen 
to be quite substantial (compared to the bare tc) which 
can be traced back to the relatively large screening k s 
(r s = 0.86). We also have t — > as k^T, £ — > 0. For very 
small r s or increasing T we can drop vv n <§; 1 in Eq. Q , 
and thus recover the Born approximation with the bare 
potential Ve- 
in addition, the higher-order terms appearing in the 
iteration of the Bethe-Salpeter equation further reduce 
the scattering (compare t to the first order V) In 
Fig- EJb) we see the significant reduction of the scattering 
length A as the density n is increased, which could be 
tested experimentally. 
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Figure 2: Plots of scattering quantities obtained from the 
t-matrix Q for r s = 0.86 (GaAs) and k B T/E F = 1(T 2 . (a) 
Angular dependence of the scattering length \(6). We com- 
pare A to its Born approximation given by t ~ V and to the 
bare scattering of two particles (no Fermi sea) given either by 
the exact result (tc) or by the first-order (Vc). (b) Depen- 
dence of X(ty/2) on the sheet density n (see the corresponding 
r s — mel/h 2 ^/Tm in the top axis), (c) TS-ratio R(9,5Q) of 
triplets/singlets detected at a scattering angle 6, for an aper- 
ture 69 = 5°. 

Production of EPR-pairs. We now turn to the pro- 
duction of entangled electrons in the spin-singlet state. 
We consider detectors placed at 0, with a small aper- 
ture angle of 280 (Fig. QJ. We introduce the scattering 
lengths around (integrated over 250) for singlets (+) 
and triplets (-): \ S/T (0) = 2 d0'\f(0') ±/(tt- 0'\ 2 , 
A useful measure is the ratio between the number Nt/s of 
detected triplets/singlet, R(0,50) = N T /N S = 3\ T /\s- 
Next, expanding for small 50 around = it/ 2 we 
find in leading order As ~ 8\f(ir/2)\ 2 50, and At — 
(8/3)|/'(7r/2)|W, which yields 



R(tt/2,50) 



/W2) 



/(t/2) 



(6) 



Using Eq. © we find |/'//| = 0.48 at = tt/2 for 
GaAs, and note that this ratio remains of order unity 
for a wide parameter range kBT/EF = 10 _1 — 10~ 10 , 
and r s = 0.1 — 1. The Born approximation |V/V| = 
1/2(1 + r s ) ~ 0.266 is approached for k B T/E F > 10" 1 . 
Therefore, this setup indeed allows the production of 
dominantly singlets (EPR pairs) at the detectors 1 and 
2 provided that the aperture angle is sufficiently small. 
For example, i?(90°,5°) = 0.2% or #(85°, 5°) = 0.7%; 
see Fig. Efc), 

To estimate the singlet current for a given input cur- 
rent I we assume that the incident electrons occupy the 
lowest transverse mode in the QPC, giving plane waves 
of transverse width w (typically w ~ 100 nm). The prob- 
ability for the singlets to be scattered into the detec- 
tors is P s = (l/A)\s/w = 0.06% with A s = 0.24 nm 
for 50 = 5°. It is advantageous to inject simultane- 
ously the two electrons from the reservoirs (e.g., by 
opening both QPCs simultaneously) (2j|. Then the 
singlet current is given by Is = Pgl = 0.6 pA for 
1=1 nA. The total scattering length for unpolarized 
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electrons is A to t = (1/4) f^ 2 /2 d9[X s (6) + 3A T (0)] = 3.4 
nm in GaAs (compared to 11 nm in Born approxima- 
tion). This is consistent with the ladder approximation 
(Atot&F = 0.53), and yields the total scattering probabil- 
ity Ptot = Atot/w = 3%. 

An important requirement is that only correlated elec- 
trons (singlets) are counted at the detectors 1 and 2, i.e. 
electrons which have scattered off each other at 9 = 7r/2, 
whereas we need to avoid the counting of uncorrelated 
electrons which are accidentally scattered into the de- 
tectors due e.g. to impurity scattering. This require- 
ment could be fulfilled e.g. by coincidence measurements 
and/or with the help of an ac modulation applied to 
each reservoir with different frequencies uji and u)2- This 
would enable a frequency selection of the electrons, since 
only the electrons which have interacted are modulated 
by the two frequencies uj\ ± 102 ■ 

Besides the observation of the scattering length and 
its density and angular dependence, our proposal could 
be further tested by adding a beam-splitter to probe 
the singlet-state via noise measurement 0- E3|- by per- 
forming tests of Bell's inequality K k . lift with spin filters 
1 1 2L 1 1 . or by using p-i-n junctions |22j to transform sin- 
glets into entangled photon pairs. An alternative test 
requires spin-filters, obtained eg. by tuning the QPCs 
into the spin-filtering regime II J . or by replacing them 
by spin-filtering quantum dots |l2j • Then with increasing 
spin polarization V the probability of incoming singlets 
ps = (1 — "P 2 )/2 is suppressed, and the singlet current at 
6 = 7r/2 vanishes at full polarization V = 1. 

We now comment on the Kohn-Luttinger instability 
[l^j . The crossed diagram — which in 2D can lead to 
an instability only for excitations with q > 2kp 23] — 
does not lead to a strong renormalization of the scatter- 
ing vertex, as the associated temperature is infinitesimal, 
k B T/E F ~ exp(-10 3 ) Q. This is larger than the value 
~ exp(— 10 5 ) found in 3D [3, despite the asymptotic 
decay for v n ~ n~ 2 being slower than in 3D, vi ~ e~ l 
(the decay is polynomial because the 2D potential is non- 
analytic). The crossed diagram, given by the polarization 
propagator, has asymptotics ~ n -3 / 2 , instead of ~ l~ 4 in 
3D. Finally, we also checked that the repulsive electron- 
electron interaction is not appreciably affected by polar 
or acoustic phonons [U 

In conclusion, we have shown that the scattering pro- 
cess of two electrons can be used to isolate and select 
the entangled spin-singlets at the scattering angle 7r/2. 
The process survives in the presence of an interacting 
Fermi sea which, however, reduces the scattering ampli- 
tude substantially. Nevertheless, according to our anal- 
ysis the current of entangled electrons is within experi- 
mental reach for realistic two-dimensional electron gases. 
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